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1 Introduction 



Quantum Chromodynamics (QCD) has been accepted as the most successful theory to 
describe the strong interaction. However, the low energy property, such as the property 
of baryon, is hard to be described with the precise QCD. Lattice QCD is then devel- 
oped to calculate the property of baryons from the full QCD lagrangian (see for example 
Ref. PQ). Meanwhile, on the side of continuous field theoretical method, there are lots 
of phenomenological models to try to give the picture of the baryons. Soliton models, 
including the two most successful models, Skyrme model(or chiral soliton model, see for 
example Ref. [2]) and the chiral quark soliton model(or Nanbu-Jona-Lasinio (NJL) Model, 
see for example Ref. [3 ), are one kind of these models to give a vivid picture of baryons. 
However, the Skyrme model and the chiral quark soliton model give two contrary pictures. 
In Skyrme model, baryons are described as a soliton of mesons with no valence quark de- 
grees of freedom. In contrast to Skyrme model, baryons are described as valence quarks 
bounded in the meson soliton background in the chiral quark soliton model. However, 
only local interactions are taken into account in NJL model. 

To search for the solid QCD foundation of the soliton models and implement the nature 
of non-locality of the quarks in mesons, the global color symmetry model (GCM) has been 
developed since the middle of 1980'sjU El EH 13 IE]- In the GCM, the gauge symmetry was 
discarded and then QCD is reduced to a finite-range current- current interaction theory. 
One can get a quark-meson interaction model or a quark-diquark interaction model after 
bosonization[S] of the current- current interaction theory. Since there are many kinds 
of meson fields to interact with quarks and it is hard to deal with all of them, one 
usually takes the low energy effective mode, the Goldstone mode pion and sigma, and 
discards the others. Then light baryons are modelled as solitons in the chiral meson 
background with Goldstone mesons PI ITU1 ITT] . With the GCM soliton, one can also discuss 
the dynamical confinement and, in particular, the effect of the non-locality. Furthermore, 
the dependence of baryon properties on hadron matter medium has also been studied in 
the GCMjUl H3J HH HSl CSl HZl HBl However, in the practical calculations, one usually 
takes the chiral circle constraint, i.e., solves the coupled equations of quark and cx-meson 
field explicitly, but discards the pion fields. Then in this paper, we extend the GCM soliton 
to beyond the chiral circle constraint, i.e., deal with the scalar field a and the pseudoscalar 
field 7r simultaneously. We can recognize that there is much difference between the results 
with and beyond the chiral circle constraint. 

Below, we will give a concise description of the GCM model and of the soliton in 
GCM model in section II (an extensive review of GCM model can be found in series 
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of papers [T§1 I2U]). In section III, we describe the method of the numerical calculation 
and give our numerical results and compare the results with those with the chiral circle 
constraint. Finally, we summarize our work and give a brief remark in section IV. 



2 Brief Review of GCM 

The global color symmetry model (GCM) is constructed by the the generation functional 
in Euclidean space asjUGI] 

Vq{x)Vq{y) exp[-J d A xq{x)fiq{x) -^JJ d A xd i yq{x)^ pL —q{x)D{x-y)q{y)^ tl —q{y)} . 

(1) 

This functional is invariant under global color SU(3) transformation rather than the 
gauge color SU(3) transformation. D(x — y) is the effective gluon propagator, which is a 
parameterized function to model the low energy property and dynamics such as those of 
hadrons. It has been shown that the infrared enhancement of D(x — y) is closely related 
to the chiral symmetry breaking (see for example Refs. |21[ I22j). Such a character should 
be embedded in the model. 

It is convenient to calculate the Grassmann integration after bosonization[4j. Firstly, 
one can reorder the quark fields through the Fierz transformation and get 
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Z = / VqDqexp 



q 
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q{x)-^-q(y) D (x - y)q{y)^q{%) 



(2) 



where M e = K a ® C h ® F c , with {K a } = {I,i l5 , j-^, ^TstJ, {C b } = {^h 7^}, 
{F c } = {-^, are the direct product of spin, color SU(3) and flavor SU(2) matrices 
resulting from Fierz transformation. Then, one introduces bi-local Bose fields B e (x,y) 
which transform like q(x)^j-q(y). The bilocal fields have the quantum numbers as those 
of mesons, so the fluctuations on vacuum can be identified as mesons. One obtains then 



Z = J VB e (x,y)e~ 



■S[B e (x,y)] 



with an action given as 



S [B e (x,y)]= Jd 4 xd 4 yq(x) 



1 -dd(x-y) + —B e {x,y) 



q{y) + / d 4 xd A y 



(3) 



B%x,y)B e (y,x) 
2gW{x-y) 

(4) 



Introducing the quark chemical potential through the transformation 



q(x) -> e tlX4 q(x) 



(5) 
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Sfr,B>] = -Tr\nG-^B e ] + / dWy ^'gf^f , (6) 



(7) 



and integrating the quark fields, one gets the action 

B e (x,y)B (y,x) 
2gW{x-y) 

where the inverse of the quark propagator can be written as 

x,y) = e flX4 G- 1 {x,y)e- fiy4 

= (7 • d - ^ 4 )5(x -y) + e» X4 ^B 9 (x, y)e~^ 4 

Generally, the bilocal fields can be written as 

B d (x, y) = B°(x, y) + j: rf (z, y)rf(^) , (8) 

where the first term is the translation invariant vacuum configuration of the bilocal fields. 
The second term is the fluctuation over the vacuum which can be identified as effective 
meson fields. Y 6 is the form factor of the meson fields, which we take as the form of B® 
with appropriate Lorentz structure since the excitation of the inner degrees of freedom 
needs much energy and can be frozen in low energy range [Hj. There are many meson 
fields, including scalar, pseudoscalar, vector, axial vector and tenser field. In the lowest 
order, we can only consider the Goldstone mode, 4> d = {a, tt}, which was thought of as 
the most important low energy degrees of freedom. 

The vacuum configuration can be determined by the saddle point condition jj=p = 
0, which induces a translation invariant quark self-energy S(x — y) — ^-B^x^y) = 
^j-Bq(x — y). The equation of S(x — y) in momentum space is, in fact, a truncated 
Dyson-Schwinger equation, which reads 



E(p) = i 1 ■ p[A(p 2 ) - 1] + B(p 
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g2 J #^- 9) Y 7 ^ 7 .gls(g) 7 4' (9) 

where D(p) is the Fourier transformation of D(x). 

It is evident that the quark meson coupling is in Yukawa form while the meson self- 
interaction is nontrivial. The dynamical chiral symmetry breaking is due to B(p 2 ) 7^ 
0, which causes quark a dynamical mass B(p 2 )/A(p 2 ). One can model the low energy 
property through a certain form of A(p 2 ) and B(p 2 ) determined by solving Eq.@ with an 
effective gluon propagator or derived from some other models such as instanton model, 
or lattice stimulation. 

In the chiral limit, the bi- local fields can be written as 

— [B 9 (x,y) - B°(x,y)} = fp x (iJ) e *»^W//-, (10) 
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where r = x — y and R = can be considered as the relative coordinate and centroid 
coordinate of a quark and a anti-quark in a meson, respectively. 

Expanding the action to leading order in derivatives of o = x cos [0//tt] and tt = 
0X s i n [0//7r] at the vacuum a = f n , n = 0, one gets 

S[a, 77] - S[U, 0] = / d'R{^(d,a) 2 + l -{d^) 2 + U{ X 2 {R))}, (11) 
where £7(x 2 ) is the effective potential of mesons with x 2 = cr 2 + 7? 2 and can be given as 



u(x 2 



-12 



d A q 



In 



g 2 A 2 (g 2 )+E 2 (g 2 )( X /A) 
g 2 A 2 (g 2 ) + B 2 (q 2 ) 



B 2 {q 2 )[{ X /U) 2 -\ 
q 2 A 2 (q 2 ) + B 2 (q 2 ) 



and is the pion decay constant given by 



U = 12 



B 2 A 2 



k 2[ 



g 2 [(5') 2 + 55"] + 55' 



B< 



(12) 



(13) 



[ g 2 A 2 + S 2]2 g 2^2 

It should be mentioned the above expression of f n is only correct when the derivatives 
of A(p 2 ) are neglected. However, this approximation can be fairly right according to the 
characteristic of the effective gluon propagator in the present calculation. The full formula 
to evaluate the f n can be found in [Zj. 

The meson masses can be obtained by differentiating the potential U twice with 
respect to a and tt. It can be found that the pion mass is zero and the a mass is finite at 
the chiral limit. 

Originally, baryons are regarded as solitons with bag constant [3] 

B = U(0) - U(U) , 

where U(0) is the effective potential of mesons in the soliton with assumption a = 0, if = 
and U(f n ) the one of vacuum with a = f n , tt = 0. Even though some baryon properties 
and their dependence on baryon density have been obtained in such a model (see for 
example Refs.jH^]), the advantage of the quark-quark interaction in the GCM has not 
yet been represented explicitly. Then a much more elaborate soliton model was developed 
in 1990 's [1(3]. As in Ref.^0], the meson fields in the soliton are written explicitly as 

M 6 



■[B"(x,y)-B°(x,y)] = B(x-y) a( 



x 



+ ?75^ ■ 7T( 



,x 



(14) 



2 y 1J/ uv v L 2 v 2 

After introducing a chemical potential of quarks as the above, one expresses the action- 
difference between the vacuum state and a system with nonzero chemical potential as 

S[fx, a, tt] = Tr In G- 1 [/x = 0, a, tt] - Tr lnG _1 [/i, a, tt]- J d 4 x[^(d x a) 2 + ^vr) 2 + U( X 2 )} ■ 

(15) 
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where \ 2 — 0-2 + The inverse of quark propagator at — is given as 

G- V = 0; x, y) = t • dj(x -y) + j-B(x - y) H^-) + , (16) 

and that at finite /x can be written as 

G-\ii;x,y) =e^G- 1 { f i = 0;x,y)e-^. (17) 

In the vacuum, the fields a is f n and fr is zero determined by the saddle point condition 
of S[a, 7r]. When there are valence quarks, the meson fields will vary with the spatial 
coordinate to show the shape of the bound system, the hadron. 

One can define an effective potential through the Legendre transformation. With the 
static meson field, the effective potential is proportional to the energy functional. At the 
Hartree level as in Ref. the energy functional of the static meson fields at a set of fixed 
quark occupation number n is 

E[n, a, 7r] = E q [n, a, tt) + E m [n, a, n], (18) 

where E q is the valence quark's contribution, 



E g {n,a, 7?] 



d 4 x 



Tr In - rrlnG _1 [0, x] -fin, (19) 



and E m the meson field's contribution 



E m [n,a, 7?] = J d 3 x 



^(Va) 2 + i(W) 2 + f/( x ) 



(20) 



The quark chemical potential \i is a functional of the meson fields a and ir and the quark 
occupation number n, which can be deduced from the relation n = dTr In G -1 [/a; a, if]/dfi. 

It is convenient to calculate the above trace of the quark propagator in Eq.(fT§)) at 
certain bases, and one can get the bases through the spectral decomposition of the quark 
propagator. With the static meson fields, G~ 1 (x,y) is time translation invariant and 
allows stationary eigenvectors of the form Uj(x)e luJX4 , which satisfies 

d 3 yG~ 1 (u; x, y)Uj{y) = z7 4 A i (w)uj(£) , (21) 

where G _1 (c<j; x, y) is the Fourier transformation of G~ 1 { y). The eigenvalues 

Xj(uj) can be written as the form Xj(u) = u — iej(u), where €j is the quark eigen-energy 
of the state j. Calculating the functional in Eq.(|19|). one can know that the total energy 
is the sum of the positive eigen-energy €j, which satisfies \(sj) = 0, up to the chemical 
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potential These positive energy states satisfy the Dirac-like equation, which can be 
given in momentum space as 

[i 7 .M(p)+5(p)] Ui (p) + j- y 7^5(^)[ff(p-fc) + i7 5 f-7f(p-fc)K-(fc) = 0. (22) 

where a = a — f n . Because of the interaction, there is a renormalization constant for the 
quark field, which can be derived from the residue at the pole of Green's function and 
given as[H] 

Z, = - / d* P d\uM dG ~ l f/: M uM . (23) 
The total soliton energy functional of the soliton is then 



3 x 



\(yo) 2 + \{vtf + u{ X ) 



(24) 

i=i LZ z J 

The static meson fields are the configuration which minimizes the energy functional, i.e., 
they are the solutions of the equations SE[a, it]/5<j = and SE[a, 7f]/<5-7f = 0. After some 
derivation, one can write these equations explicitly as 

-VV(f) + ^ + Q CT (r) = 0, (25) 

-V 2 7f(^ + -^L + Q*(f) = 0, (26) 

where the Q a and are the source terms contributed from the valence quarks, and can 
be written as 



Q*(R) = E tV / dVyu^Bi-^, x - y) ll5 rS(^^ - R) Uj {y) . (28) 

j = l J-K^j J 1 



The above Eqs. ()22|) . (|24T - (|28|) can be solved self-consistently with numerical technique. 
When one wants to solve these equations with the chiral circle constraint, one needs only 
to discard the ir degrees of freedom or set it as zero to solve the other equations. And it 
has been carried out in Refs. [01 EI] . In this work, we intend to solve the equations beyond 
the chiral circle constraint. 

3 Numerical Calculation and Results 

To solve the soliton equations, one needs at first the effective gluon propagator D(x — y) 
to get the scalar functions A(p 2 ) and B(p 2 ) and fix the quark propagator. In the present 
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work, we take the effective gluon propagator in the form of the Munczek-Nominovsky 
modeljUl 

g 2 D( q ) = (27rr^5\q). (29) 
Then the Dyson-Schwinger equation (Eq.©) yields solutions 

f 2 v 2 < ^ 

MP 2 ) = , f^r^ o ~ V (30) 



I , p 2 > \ . 

This long range interaction has been reported to have the dynamical confinement property 
|24j . Although the effective gluon propagator can take a more realistic form, we would 
show that this simple model can give a good data to compare with the experiment. 

From Eqs. (|22|) . (|24 |) -([28 |) . one knows that they are coupled Dirac-like equation in 
momentum-space and Klein-Gorden-like equations in coordinate-space. As mentioned 
above, this set of coupled equations has been solved with the chiral circle constraint, i.e., 
setting 7? = 0, in Refs.jHl [TUJ . In this work, we solve the equations with and beyond the 
chiral circle constraint with the same methods as that in Refs.|9~l ITUj. 

For the equation of quark field, we solve it in momentum space. For a nucleon, 
we consider only the states with orbital angular momentum L = 0, and express the 
wavefunction of a valence quark as 

uM = f J iiP \ , ) ■ (32) 

With the numerical integration and differentiating technique, we transform the Dirac-like 
equation to an algebraic eigenvalue equation H mn (ej)X n = 0, where X 2n = fj{nh) and 
X2n+i = 9j{ n h) with h being the integration step of the momentum, ej is the eigenvalue. 
If there exists a nontrivial solution of X, det(H) should be zero. We search for the 
eigenvalue ej from zero to a definite value until det(H(ej)) = and get the eigenvalue e 
and eigenvector X with the singular vector decomposition method(SVD)[25j. 

For the equations of the meson fields, we solve them in coordinate space with Newton's 
functional iteration method. Since the left hand side of Eqs. ()25|) and (|26[) is a functional 
of a and ff, we can refer it to F[a, 7?]. The meson fields (cr s ,7? s ) are the solutions of 
F[a s , 7? s ] = 0. To solve these equations, we expand the F[a, tt] at an initial value a = a , 
7r = 7? o and obtain 



F[<7,7?] = F[a ,7r ] + j d 



6a 



5a(x) 

<T=cro,7?=7?o 
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5F[a, tt] 
5tt 



5n(x) + • • • . 



(33) 



Neglecting the higher order and setting F[a, tt] as zero, we gain an iteration equation of 
the meson fields, 



where u n {x) and 7r n (x) are the solutions after n-times iteration. It should be mentioned 
that, with the chiral circle constraint, i.e., 7? = 0, the coupled equations can be solved with 
less effort. However, when one solves the equations beyond the chiral circle constraint, 
there are four fields: the scalar field a and the triplet tt fields. It is then very difficult to 
solve these equations. To the convenience of numerical calculation, we take the Hedgehog 
form solution 



which is one of the solution of the coupled equations. 

In the calculation, we take at first a set of trial meson fields to solve the Dirac-like 
equation and get the eigen-energy and wave-function of the quark and the renormalization 
constant. With the source terms constructed from the quark field, we solve the meson 
field equations to get the new meson fields. Inserting the new meson fields in the quark 
equation and repeating the above steps until the meson fields and quark wave-function 
converge to a desired precision, we obtain the final solutions. 

One can recognize that f n is in proportion to r] from Eq. (jl3j) when one takes the 
effective propagator in the form of delta function in Eq. (j29|) and gets the function A and 
B as Eq.()30|). One can then infer that the Eqs.()22|). ()25|) and ()26|) can be scaled with the 
dimensional constant 77. All the results can thus be scaled to dimensionless data. 

After solving the coupled equations, we obtain the eigen-energy and wave-function of 
the quark and the meson fields (the obtained results of the rescaled quark field and meson 
fields are illustrated in Fig. Q]), and further the potential energy E p and kinetic energy Ef. 
of the meson fields. The energy of the soliton can then be determined by E = 2>Ej+E p +Ek. 
Furthermore, We obtain the mass of the soliton with two methods. One is the naive 
center mass reduction, M s = \/E 2 — 3 < p 2 >, where < p 2 > is the expectation value of 
the square of the quark momentum. The other is the recoil correction [2 7j. 

In order to get a more impressive picture of the soliton, we adopt explicitly the results 
with 77 = 1.04 GeV, which is fixed by fitting the experiment data = 93 MeV^O]- The 




a n (x) + 5a n (x) 

TT n (x) + 5ll n (x) 



(34) 



a(R) = a{R) , 



3). 



(35) 
(36) 
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Figure 1: Rescaled quark and meson fields in the soliton obtained in the numerical 
calculation. The thinner and the thicker curves correspond to the results with and beyond 
chiral circle constraint, respectively. The left top, bottom panel shows the upper, lower 
component of the quark field accordingly. The right top, bottom panel displays the a, n 
field which is scaled with the pion coupling constant f n , respectively. 

obtained properties of the soliton are listed in Table 1. The obtained quark wave-function 
and the meson fields are illustrated in Fig El 

From Table Q and the Fig |21 one can recognize evidently that, besides the quite 
extended distribution, the pion field provides a strong attractive interaction and reduces 
the energy of the valence quark from 339 MeV to 129 MeV and further the whole energy 
of the soliton from 1430 MeV to 1124 MeV. From Fig one can see that the a field 
beyond the chiral circle constraint is a little weaker than that with the constraint. Then, 
the interaction between quarks with a field is weaker. However, one can get a strong 
7r field beyond the chiral circle constraint, which provide a strong attractive interaction 
between quarks. Thus, one can get a tighter bound and lighter soliton beyond the chiral 
circle constraint. 
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Table 1: Calculated mass and root mean square radius of the soliton with and beyond 
the chiral circle constraint (with the strength parameter rj = 1.04 in Eq.(|29|) ) 





with 


beyond 


energy of valence quark e (MeV) 


339 


129 


meson potential energy(MeV) 


218 


109 


meson kinetic energy(MeV) 


267 


697 


energy of soliton (MeV) 


1502 


1193 


mass (naive center of mass correction) (MeV) 


1430 


1124 


mass (recoil correction) (MeV) 


1020 


916 


radius (naive center of mass correction) (fm) 


0.71 


0.64 


radius (recoil correction) (fm) 


0.59 


0.56 



In the calculation with chiral circle constraint, we have a state with good quantum 
numbers, such as angular momentum and isospin. Then we can compare the result with 
the experiment data of a nucleon. The results listed in Table 1 show that the calculated 
result does not agree with experimental data well. In the case of beyond the chiral circle 
constraint, our calculation can not give a state with good quantum numbers because of 
the inclusion of the pion field and the Hedgehog approximation (The Hedgehog state is 
invariant under the simultaneous Lorentz rotation and the isospin rotation. In order to 
derive a state with good quantum numbers, we should quantize the classical soliton). The 
state we obtained is then a mixture of a nucleon and a delta. The result listed in Table 1 
indicates that the presently obtained mass of the soliton with the naive correction on the 
center of mass is comparable to the experimental data. However, in the case of the recoil 
correction on the center of mass, the obtained mass is smaller than the experimental data. 
Such a result is not strange since we have not quantized the soliton to a physical state. 
Meanwhile, it is promising to see that the soliton mass is close to the right. When we 
quantize the soliton, the nucleon mass is about tens of MeV heavier than the soliton mass. 
Then the mass of a nucleon is nearly equal to the experiment data. The Table shows also 
that the radius of the soliton is smaller than the nucleon charge radius 0.87 fm observed 
in experiment. It is quite natural because we have taken the effective gluon propagator 
with infinite long range interaction. The realistic interaction is surely not in infinite long 
range. If we take a more realistic form with long range decreasing behavior, we expect 
that a better fit to the experiment data can be obtained. 
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Figure 2: Calculated quark and meson fields in the soliton with and beyond chiral circle 
constraint in the case that the energy scale 7] = 1.04 GeV. The thinner, thicker curves 
correspond to the results with and beyond chiral circle constraint, respectively. The left 
top (bottom) panel shows the upper (lower) component of the quark field. The right top 
(bottom) panel displays the a (ir) field which is scaled with the pion coupling constant 

u 

4 Summary and Remarks 

In this paper, we studied the GCM soliton with and beyond the chiral circle constraint. 
The calculation shows that, when the pion field is taken into account in the scheme 
beyond the chiral circle constraint, the calculated mass of the GCM soliton is comparable 
with experimental data. The calculation indicates that the pion field provides a strong 
attraction, which binds the quarks in the soliton more tightly and decreases the energy 
of the soliton to a more realistic value. 

In this calculation, we take the form of the effective gluon propagator as a delta 
function in momentum space. In such a case, the interaction between quarks is in infinite 
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long range with a constant strength. The calculated mean square root of the radius of the 
soliton is then smaller than the experimental data. Since the realistic form of the effective 
gluon propagator should include the infrared enhancement, which is closely related to the 
chiral symmetry breaking, and the ultraviolet property of asymptotic freedom, which has 
a decrease on the interaction strength, more precise calculation with a realistic effective 
gluon propagator is needed to give a good description of the experimental data of nucleon. 
Besides, the soliton we obtained with the pion field being included explicitly does not 
possesses the good quantum numbers of a nucleon. Then, a quantization on the classical 
soliton and calculations of more physical observables are necessary. The relevant works 
are under progress. 
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